Enhanced ferroelectric polarization by induced Dy spin-order in multiferroic DyMnOs 
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Neutron powder diffraction and single crystal x-ray resonant magnetic scattering measurements 
suggest that Dy plays an active role in enhancing the ferroelectric polarization in multiferroic 
DyMnOa above T^^= 6.5 K. We observe the evolution of an incommensurate ordering of Dy mo- 
ments with the same periodicity as the Mn spiral ordering. It closely tracks the evolution of the 
ferroelectric polarization which reaches a maximum value of 0.2 ^iC/rn^. Below r°^, where Dy 
spins order commensurately, the polarization decreases to values similar for those of TbMnOs. 

PACS numbers: 61.12.Ld, 61.10.-i, 75.30.Kz, 75.47.Lx, 75.80.-|-q 
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The strong coupling between ferroelectricity and mag- 
netism in modern multiferroics has offered a new 
paradigm of magnetoelectric materials. It has stim- 
ulated a search for new multiferroics and highlighted 
the need for a deeper understanding of the physics 
behind these exciting materials. Although the re- 
quirements of magnetism and ferroelectricity are chem- 
ically incompatible in these new multiferroics spin 
frustration'2] leads to complex magnetic arrangements 
that can break inversion symmetryUQ. The strong cou- 
pling provided by frustration in the perovskite mangan- 
ites i?Mn03 leads to a number of magneto-electric phe- 
nomena. For example, compounds with i?=Tb and Dy 
exhibit flops of the direction of the spontaneous polariza- 
tion (Ps) with applied field (H) while ferroelectricity is 
observed only under a magnetic field in R—Gd^^^J^. 
Of these perovskites, i?=Dy exhibits the largest value 
of Ps=0.2 ^C/m? three times larger than that for 
R =Tb) and a giant magnetocapacitance effect^]. 

In these materials, a phenomenological treatment of 
the coupling of a uniform electric polarization P to an 
inhomogeneous magnetization M leads to a term lin- 
ear in the gradient VM, the so-called Lifshitz invari- 
ant, that is allowed only in systems with broken inver- 
sion symmetry!^. This model is consistent with neutron 
diffraction experiments on TbMnOa that showed that a 
spiral arrangement of Mn-spins within the 6c-plane de- 
velops at T; = 28 K, coinciding with the onset of PsQ. 
Although in this picture the contribution to Pg from the 
magnetic ordering of i?-ions is ignored, their role is under- 
scored as a source of magnetic anisotropy that is required 
to predict the correct direction of Pg under an applied 
magnetic field 

In this Letter we present measurements on DyMnOa 
suggesting that Dy plays an active role in enhancing the 
ferroelectric polarization in multiferroic DyMnOa above 
T^^= 6.5 K. We have combined neutron powder diffrac- 
tion and single crystal x-ray resonant magnetic scatter- 
ing to investigate the evolution of magnetism in DyMnOa 
with temperature. We find that the Dy moments order 



in an incommensurate (ICM) structure with the same pe- 
riodicity as the Mn moments below a temperature T^^— 
15 K. The transition from the commensurate (CM) Dy 
ordering below T^^ to the ICM state is associated with 
an enhancement of Pg just above T^^ to a value ap- 
proximately twice that found for the maximum Pg in 
TbMnOa. The CM-ICM transition shows a large hystere- 
sis in which Pg and the intensity of magnetic reflections 
arising from the ICM Dy magnetic ordering exhibit a 
similar behavior. Our work suggests a magneto-strictive 
coupling between Mn and Dy spins giving rise to an en- 
hanced Ps above T^^ compared to other perovskite mul- 
tiferroics. This mechanism would add to the ferroelectric 
polarization that arises from the Mn spin-spiral. 

In the multiferroic perovskite manganites it is found 
that for i?=Tb below TAr~41 K, Mn-spins order first 
along the 6— axis in an ICM sinusoidal arrangement with 
propagation vector t^^" = (0 0.28. ..0.29 0)|1 [l^. A 
similar behavior with r*^" = (0 0.36. ..0.385 0) is ex- 
pected from x-ray measurements for DyMnOsj^]. Sec- 
ond harmonic lattice reflections {q — 2t) associated with 
the magnetic ordering have been observed to arise from 
a coupling of the ICM magnetic ordering to the lattice 
via a quadratic magneto-elastic coupling pdL Il2|. Below 
T/=28 K for TbMnOs, an additional component of the 
Mn magnetic moment along the c-axis, in phase quadra- 
ture with the component along b, gives rise to a spiral 
magnetic ordering and breaks inversion symmetry, lead- 
ing to the observation of Ps along the c— axis. In this 
regime Mn spins induce an ICM ordering of Tb-spins 
with the same propagation vector]^ . A similar transition 
occurs for i?=Dy at r; = 18K0 as illustrated in Fig. 1(a). 
Below T§ <10 K, Tb and Dy magnetic moments order 
separately with propagation vectors T"^* = (0 0.42 0) and 

TDy={o i o)E 14]. 

To investigate the evolution of the magnetic order- 
ing of DyMnOa with neutron powder diffraction (NPD) 
we used a 0.65g polycrystalline ^^^DyMnOa sample pre- 
pared from a mixture of isotope enriched ^^^Dy203 ox- 
ide (94.4% enrichment) and Mn20a, using standard solid 
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FIG. 1: (Color online) Temperature dependences of (a) elec- 
tric polarization along the c-axis for DyMnOs and TbMnOs 
(from Goto et al.^), (b) integrated (neutron) intensities of 
magnetic A-type and (0 | 1) reflections measured on heating 
(I ™'Mn3+ shows the intensity corresponding to a 4 /ib Mn"^"'' 
magnetic moment for this reflection; shaded area corresponds 
to the magnetic intensity on top of the Mn^"*" saturated mo- 
ment), (c) temperature dependent integrated (synchrotron) 
intensities of selected Dy commensurate (heating) and incom- 
mensurate (heating and cooling) superlattice reflections. The 
region around T^^'is shown in the insert. Here, AI highlights 
the difference between intensities on cooling and heating at 
7 K. 



state synthesis methods. Isotopic ^^^Dy was chosen for 
its smaller neutron absorption cross section (ua) com- 
pared to that of natural Dy|0|. NPD data were mea- 
sured from this sample between 2-300 K on the GEM 
diffractometer at the ISIS-facility, Rutherford- Appleton 
Laboratory. The data were analyzed with the FuUProf 
refinement package jT^. Synchrotron x-ray diffraction 
measurements on a DyMnOs single crystal were con- 
ducted at the 7 T multipole wiggler beamline MAGS, op- 
erated by the Hahn-Meitner-Institut at the synchrotron 
source BESSY in Berlin. Details on the crystal growth, 
beamline and experimental procedure have been reported 
previously [ij, [13 • Low field (500 Oe) magnetization of 
the polycrystalline and single crystal samples was mea- 
sured in a SQUID magnetometer between 4 - 100 K. 



The data showed T, 



Dy 



N 



9 and 6 K for the polycrys- 
talline and single crystal sample, respectively. The value 
of of our single crystal is in good agreement with 
published dataQ and may suggest that the higher value 
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FIG. 2: (Color online) (a) Neutron diffraction patterns ob- 
tained from the 24°-45° bank of the GEM diffractometer at 
2, 10, 30 and 50 K. (b) fc-scans along (0 k 4) at various tem- 
peratures measured with a synchrotron x-ray energy close to 
the Dy-Ls absorption edge. No polarization analysis was em- 
ployed here in order to allow for simultaneous observation 
of structural and magnetic superstructure reflections. Inset: 
temperature dependence of the lattice modulation wave vec- 
tor obtained from a superstructure reflection of type (0 fc 3). 



obtained for the powder sample may arise from a small 
non-stoichiometry. This assumption is corroborated by 
a difference in the lattice constants between our powder 
sample and the single crystal at 300 KT4iJ. 

Turning first to the NPD measurements we found that 
^^^DyMnOa crystallized with the orthorhombically dis- 
torted perovskite structure (space group Pbnm). Cool- 
ing the sample below T/v ~ 40 K we find a series of 
magnetic satellites arising from Mn spin-ordering (see 
Fig.|2Ia)) similar to the observations of Quezel et al. [lOj 
on TbMnOs. Using their notation, these {hkl) ±r reflec- 
tions are characterized as A-type where h + k — even and 
I = odd and t — [0 t*^" 0) is an incommensurate prop- 
agation vector along the 5*-axis. An increase of satel- 
lite intensities (and t*^") with decreasing temperature 
is accompanied by the appearance of G-type satellites 
(h + k = odd and I — odd) below ~15 K, where a non zero 
value of Ps is observed as shown in Fig.ma)[^0|- Below 
T^^ the NPD data show an intensity reduction of the 
A-type satellites, which coincides with the appearance of 
a CM Dy magnetic order with propagation vector t^^ = 
(0 5 0) (see Figs.nt, |2t)0- This transition coincides 
with a sharp decrease of Pg ■ The G-type reflections are 
too weak for a quantitative comparison. Rietveld anal- 
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ysis of the NPD data between r°^< T <13 K on the 
assumption of pure Mn spin-ordering (see below) leads 
to unphysically large moments for Mn^+ (> Ajj^B /Mn) 
(Fig. ^b)). The rapid increase of the magnetic intensity 
below 13 K (shaded area in Fig. ^Jb)) and its decrease 
below suggests an additional magnetic contribution 
to the intensity of these reflections from Dy. As can be 
seen in Fig. [2Ia), the magnetic reflections are relatively 
broad and below (Fig. 3, caption) we shall argue that 
this probably arises from a distribution of the propaga- 
tion vectors. 

In order to probe directly the magnetic contribution of 
Dy to the intensities of these ICM magnetic reflections, 
we have used single crystal resonant x-ray scattering at 
an x-ray energy of 7.794 keV, slightly above the Dy-La 
absorption edge. With this resonant condition, a survey 
of reciprocal space was carried out at different tempera- 
tures. Fig.EIb) shows scans along the (0 k 4) direction as 
examples. At 16 K only the Mn induced structural mod- 
ulation q*^" is observed, producing the (0 2.77 4) Bragg 
reflection. At 8 K a pair of relatively strong additional 
(0 2.385 4) and (0 2.615 4) reflections appears, consistent 
with the magneto-elastic coupling g*^" — 2t*^". At 4.4 K 
the latter reflections have disappeared, in contrast to the 
associated peak from the Mn induced structural modu- 
lation that decreases in intensity and is shifted towards a 
q^" value of 0.81. The extinct reflections are replaced by 
significantly stronger Bragg reflections at (0 2.5 4) and 
(0 2.905 4) which correspond to the CM magnetic order- 
ing of the Dy moments and the associated incommensu- 
rate lattice modulation q^^ = 0.905 2t^^ respectively. 
This transition has been discussed recently elsewhere [l^. 

In order to identify the nature of the Bragg reflections 
associated with t^" above T^^ we employed linear po- 
larization analysis^. Figure 01 shows the dependence 
of the intensities of both the (0 2.385 4) and the related 
(0 2.77 4) reflections on the polarization analyzer conflg- 
uration (cr a' vs. a tt') at 8.5 K. The characteristic 
behavior shows that the former is of magnetic and the 
latter of structural origin. Tuning the x-ray energy to a 
value 20 eV below the Dy-La absorption edge leads to 
a reduction of the intensity of the magnetic (0 2.385 4) 
reflection measured in ct — > tt' configuration by a fac- 
tor of 40. This resonance enhancement shows that the 
(0 2.385 4) reflection is due to ordered Dy magnetic mo- 
ments and that any contribution of the ordered Mn mo- 
ments to this reflection is negligible. 

Figure Q^c) shows the temperature dependence of the 
integrated intensities of two particularly strong Bragg 
reflections, (0 1.385 2) and (0 1.5 2), related to r^" 
and T-^^, respectively. The half-integer reflection, as- 
sociated with the CM magnetic ordering of the Dy mo- 
ments, vanishes above T^'' = 6.5 K. Simultaneously, the 
intensity of the ICM reflection increases steeply on heat- 
ing. At the same temperature Ps||c increases to approxi- 
mately twice its value compared to that at temperatures 



(0 2+t"" 4) 
magnetic 




2.40 2.72 2.74 2.76 2.78 2.80 2.82 
k{K.\.\i.) 



FIG. 3; (Color online) Linear polarization analysis of (left) 
the (0 2.385 4) and (right) the (0 2.77 4) superstructure re- 
flections at 8.5 K. The intensity variation as a function of the 
polarization analyzer configuration shows that (0 2.385 4) is of 
magnetic and (0 2.77 4) of structural origin. The line widths 
of both reflections are signiflcantly enhanced with respect to 
the instrumental resolution that is marked by horizontal bars. 
Since the line width of the second order reflection (0 2.77 4) is 
roughly twice as large as that of the flrst order (0 2.385 4), the 
broadening is more likely due to an inhomogeneous distribu- 
tion of rvalues in the sample rather than to reduced coherence 
lengths. 



just below (Fig- Ea)). On further heating, the 

(0 1.385 2) Bragg intensity passes a maximum around 
8 K, decreases monotonically with a concave curvature 
- typical for an induced magnetic moment, and flnally 



vanishes above T, 



Dy 



15 K. In addition, it exhibits a 



signiflcant hysteresis (see inset in Fig. ^c)). On cooling 
from 16 K, the (0 1.385 2) Bragg intensity measured at 
7 K is about twice as large as the intensity obtained on 
heating from 4.5 to 7 K. Monitoring the count rate at 
flxed k — 1.385 while sweeping the temperature at vari- 
ous rates (0.05.. .1 K/min.) we verifled that whenever one 
does not cross , the cooling and heating curves above 
are reversable and the hysteresis is independent of 
the heating/cooling rate. This behavior of Dy induced 
ordering strongly resembles the hysteresis of the electric 
polarization shown in Fig. ^a), where a factor of 2 dif- 
ference is shown in the cooling and heating curves at 7K, 
the same difference that we flnd in the hysteresis of the 
intensity of the (0 1.385 2) reflection. Finally, the insert 
in Fig. |2Ib) shows the temperature dependence of g*^" 
measured in two successive heating and cooling cycles 
between 7.5 and 40 K. Clearly, a signiflcant hysteresis of 
about 5 K in width is observed in the temperature de- 
pendence. Neither on cooling nor on heating any sharp 
lock-in transition is observed. The edge of the heating 
curve, however, is consistent with the reported value for 
Ti = 18 K. 

Having directly established the behavior of Dy-spins 
above and below T^"^ we turn our attention to models 
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of magnetic ordering for Dy that can be obtained from 
the analysis of the NPD data. In the CM regime below 
T^^ the diffraction data indicates a cell doubling along 
b arising from antiferromagnetic Dy spin-ordering. Ri- 
etveld analysis shows that the Dy moments with abso- 
lute value TO^^ — 6.2(6)/iB/Dy He in the ab plane, being 
canted 30(±10)° from the b axis, and are stacked antifer- 
romagnetically along the c-axis. 

In the analysis of the induced magnetic ordering of 
Dy-spins at 10 K, above , we have taken solu- 
tions of the magnetic structure that are allowed by 
the T2 (81 representation of the Pbnm space group 
indicated b y t he spiral arrangements of Mn-spins in 
TbMnOala UJ- In the case where two irreproducible 
representations are coupled0|, the R moments are al- 
lowed to have components along the three principal crys- 
tallographic directions denoted as m^^ =(aa;, fy, fz)[l8j. 
Our Rietveld analysis shows that the magnetic moment 
is =(0.4(6),2.5(2),0.0(3)) /lij/Dy. The maximum 

value of the moment here is approximately twice as much 
as that found for the induced moment of Tb in TbMnOa 
(m'^^ =(1.2,O,O)^B/Tb0). In these refinements we as- 
sumed the same magnetic symmetry for Mn as that re- 
ported for TbMnOslil. 

The CM structure of Dy-spins below T^^ preserves the 
mirror plane perpendicular to the c axis and thus obvi- 
ously does not allow for a contribution to along the 
c— axis. Polarization measurements show that this spin 
ordering also does not contribute to along the other 
principal crystallographic directions (P||a~P||6~ 0)0. 
For the region just above our data is consistent with 
a sinusoidal modulation of the Dy-spins along the 6— axis. 
The absence of a c— axis component would exclude the 
presence of a Dy 6c-spiral as found for the case of Mn 
spin-ordering in TbMnOa. Thus, on the basis of symme- 
try the induced Dy-spin order alone is unlikely to result in 
a ferroelectric polarization along the c— axis above . 

Compared to i?=Tb or (Eu,Y), the behavior of 
DyMnOa is unusual, as Ps just above T^^ is three times 
larger for comparable conditions and shows a large hys- 
teresis. Indeed Goto et al.^ find almost a 50% difference 
in Ps at 7 K obtained after cooling to 2 and 7 K, a be- 
havior that is similar to the hysteresis of the magnetic 
intensity from the induced Dy spin-order (Fig.Q^c)). Be- 
low , Ps decreases to a value of 0.06/iC/m^ which 
is similar to values found for non-magnetic i?-ions and 
i?=Tb at 2K. This would indicate that in this regime fer- 
roelectricity arises from a similar coupling in all of these 
materials. 

Our measurements indicate that there is an additional 
(or alternate) mechanism by which Dy-spins can provide 
a significant enhancement of Ps above the value expected 
from a Mn spiral alone. Although our measurements can 
not provide a detailed microscopic model of the lattice 
distortions that arise from the magnetic ordering of Mn 
and Dy, we can exclude a mechanism based on a be Dy 



spin-spiral. The temperature hysteresis of Pg and the in- 
duced Dy magnetic reflection occur over the same region, 
suggesting a close coupling between these two quantities. 
The lattice distortion that arises from the Dy induced or- 
dering is significantly larger in magnitude and different in 
nature than those for i?=Tb. The hysteretic behavior of 
Ps and the intensity of the ICM Dy magnetic reflections 
suggests that the enhancement of P, above arises 
from a magneto-strictive coupling between Mn and Dy 
spins. Indeed this may be a peculiarity of the single ion 
anisotropy of Dy since it is a Kramers ions and Tb in 
contrast is not. Clearly the magneto-strictive lattice dis- 
tortions must be conform to r2 symmetry. 

The present results are important for a full descrip- 
tion of the magneto-electric coupling in perovskite man- 
ganites. While P-spins thus far have been perceived to 
be spectators to ferroelectricity, here we show that in 
DyMnOa the induced magnetic ordering of Dy coincides 
with a region of enhanced polarization compared to other 
PMnOa manganites, while below when Dy spins or- 
der commensurately Ps is sharply reduced. We argue 
that these observation suggest a peculiar coupling of the 
lattice that enhances Ps when Dy-spins are ordered in- 
commensurately above . 
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